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Abstract. We point out that ignoring nonlinear effects in finite size scaling 
may lead to errors in estimates of the critical temperature and Binder cumulants. 
We show that the order of magnitude of these effects can be estimated from 
data at relatively small volume. Using this estimate, we propose to use linear 
fits in increasingly small temperature regions as the volume is increased (rather 
than using a fixed temperature interval). The choice of the exact coefficient of 
proportionality can be optimized and reveals interesting crossing patterns among 
estimates. We show that the new procedure works very well for Dyson' hierarchical 
model. We discuss applications of the method for 3 dimensional spin models and 
finite temperature lattice gauge theories and comment on the nonlinear effects for 
existing calculations. 



PACS numbers: 11.15.-q, 11.15.Ha, 64.60.an, 64.60. De 

1. Introduction 

The study of the Binder cumulants for isolated systems of various sizes has been an 
important tool to determine the critical temperature and universal quantities of spin 
models [XJ [21 [31 HI El El El El [SI and lattice gauge theory models at finite temperature 
[TU1 [TT1 [TZl [T51 [Til [T51 [TBj . The basic idea is that the intersections among Binder 
cumulants curves for systems of different sizes obey scaling relations often referred to 
as Finite Size Scaling (FSS). 

Linearization in the scaling variables is often used to derive approximate FSS 
relations for the Binder cumulants and the critical temperature that allow, in principle, 
simple linear extrapolations to infinite volume [HE]. However, in practice, the errors 
bars on the intersections are often of the same order as the finite size variations and 
the linearity of the data is not always obvious. Examples of such graphs are figure 7 
in [3j or figures 3, 4 and 6 of [16]. For this reason, authors will sometimes prefer to 
rely on a single intersection for the two largest volumes available that has potentially 
smaller errors [51 [TB] . 

A possible source of systematic error in these infinite volume extrapolations 
are the nonlinear terms that have been neglected in the derivation of the linear 
FSS formulas. If these effects are small enough, it is difficult to disentangle them 
from the statistical errors. Typically, we are talking about a 1-2 percent effect. 
This may sound small, however, for the O(N) hierarchical sigma models where very 
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accurate calculations can be performed, the differences between the Binder cumulants 
at successive values of N are less than 2 percent for low values of N [T7] . Consequently, 
it is reasonable to set accuracy goals of less than 1 percent for the Binder cumulants, 
especially if the estimate is used to establish that a particular model belongs to a 
known class of universality. 

In this article, we discuss quantitatively the effects of the nonlinear terms in 
FSS with an example and we propose a new procedure to obtain infinite volume 
estimates that are minimally distorted by the nonlinear effects. The procedure can 
be used for spin models or lattice gauge models and (3 is a notation that can be used 
interchangeably for either the inverse temperature in Ising-like models or 2N c /g 2 in 
gauge theory. 

One purpose of the article is to provide a specific improvement to the common 
practice of using data at the same, equally spaced, values of (3 for all value of the 
linear size N (we call this the "fixed interval procedure" ) . Clear examples of the fixed 
interval procedure can be seen in figure 6 of [2J, figure 11 of [J], figure 7 of [TUj, figure 
1 of [11] and figure 1 of [16] where the fourth order cumulant is plotted versus /3 or the 
temperature. In the large volume limit, the cumulant becomes a function (that we will 
call /) of the reduced variable — /3 C ) / '/3 c )iV 1 / 1 ' , with v the usual correlation length 
exponent. This function / can be approached at sufficiently large volume by plotting 
the cumulant versus this reduced variable. When this is done, the data "collapses" 
up to some finite size corrections. It seems clear that in order to reduce the nonlinear 
effects, we should shrink the range of (3 as N increases. 

This idea seems simple, however its practical implementation is non-trivial. First, 
we need to estimate the size of the nonlinear effects. Second, we only know j3 c with 
a finite accuracy and the interval cannot be shrunk to a size much smaller than 
the uncertainties on /3 C . We need to find a compromise between two conflicting 
requirements. On one hand, the interval needs to be small enough to reduce the 
nonlinear effect. On the other hand, if the interval is too small, we may not cover 
properly the region where the linear approximation is valid which could result in 
errors for the intersections. The main technical question that we address here is: can 
we find an optimal width for the interval were both requirements can be reasonably 
satisfied? It is difficult to answer this question for the gauge or spin models discussed 
above because it is hard to extract the nonlinear effects from the statistical errors. 
In addition, we have no independent estimates of the infinite volume cumulants and 
calculations at large volume are CPU demanding. 

We will address this question in a model where the practical issues mentioned 
above can be avoided, namely Dyson hierarchical model [18] with an Ising measure. 
In this model, it possible to block spin very accurately with methods [El [20] [21] that 
have been recently reviewed in [22] . This allows us to make very accurate calculations 
at very large volume (linear size up to 10 6 sites) that cannot be reached with ordinary 
Monte Carlo methods. The free parameter of this model has been fixed in such a way 
that free Gaussian fields scale like in 3 dimensions. The technical details are given in 
|Appcndix A[ These large volume numerical calculations will be used later to check 
the accuracy of estimates made using data at much smaller volumes (linear size up to 
256 sites as in existing Monte Carlo (MC) calculations). 

Our presentation is focused on the estimation of the infinite volume critical 
temperature and the (universal) fourth order Binder cumulant. This quantity is 
defined precisely in section [2] where we introduce other notations and state the 
problem. We provide a parametrization of the nonlinear effects for the fourth order 
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cumulant in terms of two new parameters (see equation (|2.7p ). Accurate numerical 
values of these parameters are given in |Appendix A| in the case of the hierarchical 
model. We then discuss the linearization and graphs testing its validity. 

The new method that we propose here (the "shrinking interval procedure") is 
presented in section [3l It proceeds in two steps. First, we need to estimate the 
size of the nonlinear effects using data at relatively small volume. We explain that 
the important quantity is the ratio of the linear to quadratic leading amplitudes 
(called j\j j-i). We show in |Appendix B| that this quantity can be estimated with 
reasonable accuracy using a pencil and a ruler on a graph showing the Binder cumulant 
versus /3 at relatively small volume. We then propose to use a shrinking range 
A/3 = e/3 c (/i/ f2)N~ 1 / v where e is an adjustable small parameter which we shall see, 
can picked in an optimal way by localizing crossings in infinite volume extrapolations. 

In section |4] we compare the shrinking interval procedure with the fixed interval 
procedure. The calculations are done using a number of values of (3 and linear sizes that 
are typical in existing MC calculations. We emphasize that in doing these calculations, 
we did not use our prior knowledge of the critical temperature, Binder cumulant 
or amplitudes given |Appcndix A[ however we used our knowledge of the critical 
exponents v and w. We first used a fixed value of e (0.05) and found much better result 
than with the fixed interval method. We then compared infinite volume extrapolations 
for /3 C and the cumulant, for different values of e. The results are displayed in figureS] 
where we see remarkable crossings (among extrapolation curves for different maximal 
volume) at approximately the same value of e for the two quantities. We compared 
the results for this optimal value with the accurate ones and found relative errors of 
less than 10 -5 for /3 C and of less than 5 x 10 -3 for the cumulant. 

FigureUis the most important result of the paper. It shows that an optimal value 
of e can selected using the numerical results. It is crucial to realize that the crossing 
means that for the optimal value of e, results at not too large volume provide very 
accurate results. One should also appreciate the number of large volume calculations 
involved in producing the two graphs that would take a prohibitively long time for 4 
dimensional lattice gauge theory. 

In section [SJ we discuss the effects of the nonlinear terms for the calculations done 
in the literature and we conclude with possible applications of the new method. 

2. Statement of the problem and notations 

The new method proposed in this article can be used for two types of models. First, 
the spin models where (3 denotes the inverse temperature in appropriate units. In the 
literature on the 3 dimensional Ising model [9] (3 is often denoted K . Second, the 
lattice gauge models with gauge group SU(N C ) at finite temperature T = (A^a) -1 
where N T is the number of sites in the Euclidean time direction and a the lattice 
spacing. For these models (3 — 2N c /g 2 is not the inverse temperature. However, if we 
use the one-loop scaling [IT] . 

(T - T c )/T c ~ (/? - (3 c )12tt 2 /11N 2 . (2.1) 

This formula illustrates that for the finite temperature gauge model, T > T c implies 
[3 > (3 C because the ordered phase of these models corresponds to high T and high (3. 

In summary, the notation (3 can be used for both types of models with a consistent 
meaning: the ordered phase corresponds to (3 > [3 C . For this reason, we will define a 
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reduced f3 (or K) variable: 

K=(fi- C )/0 C (2.2) 

Other choices such as (j3 — f3 c )//3 lead to similar linear behavior near j3 c but have 
different nonlinear behavior. These considerations are important if we want to connect 
with other expansions [531 HI] ■ 

We now consider finite size scaling for isolated blocks in 3 dimensions. The spin 
models are defined on symmetric cubic lattice iV 3 sites. The gauge models are defined 
N T x A" 3 lattices. We use the convention N a — N for the variable number of sites 
(while N T is kept fixed) in order to have unified notations. In both cases, the system 
is not in contact with a larger system but isolated and defined with suitable boundary 
conditions. 

Under a RG transformation where the lattice spacing a £a, we have N — * N/t 
and Ui — > £ Vi u.i , where the itj are the nonlinear scaling variables which transform 
multiplicatively. We denote u K the only relevant scaling variable (we will not deal 
with external fields) and we will only consider the effect of the first irrelevant variable 
denoted u\. We assume that the exponents y K = l/v, y\ = —A/v = —lu arc known 
with good precision as it is the case when we are testing the hypothesis that a particular 
system belongs to a well studied class of universality. In addition, we assume the 
expansions 

u K = k + u^k 2 + ... (2.3) 

ui =uf } +u^ } n+ ... (2.4) 

We define the fourth Binder cumulant as 

B 4 =< m 4 > / < to 2 > 2 (2.5) 

with m the average of the order variable (spin for Ising, Polyakov's loop in the time 
direction for gauge theory) in 3 dimensions. Directly related quantities appear in the 
literature such as U = 1 — B4/3 or Q = I/-B4 etc. Due to the lack of consensus, 
we have picked a form that has limits that are easy to remember. Both numerator 
and denominator are unsubtracted averages and so this quantity can be calculated at 
large volume without running into loss of accuracy problems occurring for subtracted 
averages. In the symmetric phase, < to 4 > is dominated by 3 < m 2 > 2 and the the 
correction is suppressed by one power of the volume (the coefficient is a renormalized 
coupling constant). Consequently, if /3 < /3 C (k < 0), B4 tends to 3 in the infinite 
volume limit. On the other hand in the ordered phase, < to 4 > is dominated by 
< m 2 > 2 and the limit is 1. At fixed N, for sufficiently small (3, B4 is close to 3 and 
for sufficiently large j3, B4 is close to 1. As N increases, the transition sharpens as 
illustrated in figure [TJ Finite size scaling [TJ |Hj implies that 

S4 = /(u«JV*" s «iJV*S...) (2.6) 

Since y\ < 0, when the volume is large enough, B4 can be approximated by a function 
of kN Vk only. If we plot B4 versus k7V !/ ' s , the curves at various N approximately 
collapse. This is illustrated in figure [T] where one can see that the larger violations of 
this approximation are observed at smaller volume. The most common procedure to 
study the intersections of the various curves at fixed N is to linearize equation (|2.6[) 
and try to extrapolate the results at infinite N. In the rest of this article, we will 
consider the effects of nonlinear terms on this procedure. Our main assumption will 
be that 

Bi(J3, N) ~ S 4 (/3 C> 00) + f lK N^ v + f 2 K 2 N 2 ^ + (co + c^N^^N^ . (2.7) 
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Figure 1. B 4 versus f3 (left) and versus kN 1 ^ (right) , for N = 4, 8, 16, 32, 64 
and 128, for the Ising hierarchical model. 



Note that we have not included terms of the form n 2 N 1 /' J which should appear as 
a consequence of the presence of k 2 terms in the expansion of u K because they are 
suppressed by a factor N~ x l v compared to the term (nN 1 ^) 2 . We have also not 
included terms of order N~ 2lj and N~' Va '. 

In the linear approximation (/2 = c% = 0), we recover the standard linear FSS 
formula for the point of intersection denoted (/3*(N, N'), B%(N, N')) between the two 
curves B±{{3,N) and B±((3,N'), namely 

0*(N,N') =/3 c + /3 c (c // 1 )L(Ar,JV / ) , 

B2(N, N') = B 4 + CqM(N, N') , (2.8) 

with 

L(N, N') = (N~" - N'-")/(N a / v - AT 1 /") , 

M(N, N') = (N-"- 1 '" - N'- UJ - 1 / l, )/(N a ^ - N 1 /") . (2.9) 

These formulas have been written in a way that makes the symmetry under the 
interchange of N and N' obvious. They may look more familiar if we replace N' by 
bN and factor out the powers of AT. Since we assume that v and uj are known in good 
approximation, linear fits can used to determine the remaining unknown parameters 
in equation (|2.8|) . Graphs based on these formulas can be found in [3] (for 6 = 2), [To] 
and in figures [2] and [3] below. 

For reasons explained in the introduction, we will make model calculations for 
Dyson's hierarchical model. For this model, the local potential approximation is exact 
and the change in the local measure under block spinning can be calculated very 
accurately at very large volume and very close to (3 C [UH [20l [21] (see [22] for a recent 
review). The free parameter of this model, usually denoted c has been fixed in such 
a way that free Gaussian fields scale like in 3 dimensions. With this choice, as the 
blockspinning reduces the number of sites by a factor 2, the linear sizes are of the form 
JV = 2 n//3 . The details of the calculations and accurate estimates of the parameters 
entering in equation (|2.7[) are given in Appendix A 
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3. The shrinking interval procedure 

In this section, we present the shrinking interval procedure advocated in the 
introduction. It may be useful to state all the steps that need to be followed. First, 
we perform linear fits of as a function of (3 at fixed volume in a sufficiently small (3 
interval. We then use the linear fits at different volumes to determine the intersections 
called P*(N,N') and B%(N, N 1 ) in equation ([278]) . Finally, we select a set of pairs 
(iV, N') and perform linear fits to determine the unknown coefficients in the linear 
FSS formula equation (|2.8p . The final result is an infinite volume extrapolation for (3 C 
and £> 4 . 

We now discuss the first step. We need to specify the (3 interval (in other words, 
its center and width) for calculations of B4 at a given volume. Given a linear size 
N, we propose to use our best estimate (3 C of (3 C obtained from smaller sizes. This is 
an iterative procedure. We need to start with some reasonable estimate (for instance 
obtained with the finite interval procedure at small size) and then keep improving. We 
propose to restrict the calculation of B\ to the interval parametrized in the following 
way: 

|/?-/3 c | <e(hlh)p c N- 1 /" . (3.1) 

We do not need a very accurate value for f\ff%. We show in |Appendix B| that it is 
easy to estimate the order of magnitude of /1//2 from a graph such as figure [T] using 
a pencil and a ruler. The value of e needs to be chosen carefully. On one hand, we 
need e small enough in order to control the nonlinear effects. On the other hand, if e 
is very small, we need a correspondingly good estimate of (3 C - In addition, when e is 
too small, the intersections may be far away from the regions where we have values of 
£>4. These two effects could in principle compensate. Unfortunately, in the example 
where we have done accurate calculations, they go in opposite directions. Namely, the 
values of /? c obtained from small volume data are below the true f3 c while the small 
volume intersections are above (3 C as can easily be seen from figure [T] 

We want to figure out if the two conflicting requirements can be partially satisfied 
for some optimally chosen value of e. The errors on (3*{N, N') and B\{N, N') depend 
not only on e but also on (3 C . The choice of e is clearly a difficult optimization problem 
and it is useful to make numerical experiments. In the following section we show that 
values of e ~ 0.1 lead to results that compare very well with the results of |Appendix A| 
and that a reasonable compromise between the two requirements discussed above is 
possible. 

We should also mention that the use of linear fits to determine the intersections 
could be replaced by more sophisticated search methods. However, in calculations 
that are CPU demanding, and where we have only data for a few values of /3, linear 
fits is the most common practice. 

4. Comparing the fixed and shrinking interval procedures 

In this section, we compare FSS results using the shrinking interval procedure 
discussed above and the fixed interval procedure where we use fixed values of (3 for 
all values of N. The calculations were made for the Ising hierarchical model. For the 
fixed interval calculation, we took the same 9 values of (3: 1.176, 1.177, 1.184 for 
all possible N between 8 and 256, or in other words, 9 < n < 24 with n the number 
of blockspinnings. For the shrinking interval procedure, we started with the fixed 
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Fixed Interval 




Shrinking Interval 




Figure 2. Empirical values of (3*(N, N') obtained with the fixed interval 
procedure (left) and with the shrinking interval procedure (right) versus L(N, N') 
for 4 sets of 6 pairs of values defined in the text. The solid line is the linear fit for 
the set (8,32) (left) and the set (64, 256) (right). The dash line is the behavior 
expected from equation 1 12. 81 1 and the accurate values of |Appcndix A| 



interval data with N < 16 and an estimate of f3 c ~ 1.1772 corresponding to this data. 
We then followed the procedure described in section [3] with e = 0.05. This value of 
e was selected from the compromise of having linear fits that were not blatantly off 
the data and at the same time keep most of the intersections within the fitted range. 
We also picked 9 values of (3 within the specified range. At each N, we used all the 
possible intersections involving the 5 closest sizes to determine f3 c using equation (|2.8[) . 
This means that we made linear fits with 15 points and the results are quite robust 
under small changes in the procedure. 

Linear fits were performed for each N in order to express B4 as a linear 
function of 0. For each pair (N,N'), it is possible to determine the intersection 
(f3*(N, N'), B\{N, N')) of the corresponding lines. These empirical values are plotted 
versus the calculable values L(N, N') and M(N, N'), defined in (|2~9|) . in figures [2] and 
[3] In order to make the graphs readable and the dependence on N and N' clear, we 
have displayed 4 sets of 6 intersections denoted (N, N') in the graphs. The meaning of 
(iV, N') is that we take the intersection between N' and all the possible lower values 
starting at N. For instance, (8,32) is a short notation for the 6 intersections: 8 and 
32, 8 x 2 1 / 3 and 32, and, 8 x 2 5 / 3 and 32. 

For the fixed interval procedure, one can see that values corresponding to the 
smallest volumes, namely the (8,32) set, the behavior is approximately linear. The 
extrapolations at infinite volume for the fit with the six points of this set are 1.17848 
for P c and 2.53435 for B4, which is not too far from the accurate values but not very 
accurate either. However, if we increase the volume the linearity and the accuracy of 
the extrapolations degrade rapidly. 

On the other hand, with the shrinking interval method, the accuracy seems to 
improves with the size. In both graphs, the fit made with the (64, 256) set is hardly 
distinguishable from the accurate result. In order to give a more detailed idea of the 
effects on the extrapolated values, the numerical values of the parameters entering in 
equation (]2.8[) are shown in table [T] for various sets of data. 
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Figure 3. Empirical values of B£(N,N') versus M(N,N') obtained with the 
fixed interval procedure (left) and with the shrinking interval procedure (right) 
for 4 sets of 6 pairs of values defined in the text. The solid line is the linear fit for 
the set (8,32) (left) and the set (64, 256) (right). The dash line is the behavior 
expected from equation 1 12, 81 1 and the accurate values of section |Appendix A| 



n 


0c 


PcCo/fl 


B 4 


CO 


15 


1.1788071 


3.74 


2.5245 


- 0.761 


16 


1.1789292 


3.54 


2.51486 


- 0.715 


17 


1.1789933 


3.40 


2.50724 


- 0.673 


18 


1.1790389 


3.28 


2.49907 


- 0.621 


19 


1.1790832 


2.99 


2.48794 


- 0.534 


20 


1.1790042 


3.51 


2.50940 


- 0.697 


21 


1.1790198 


3.35 


2.50336 


- 0.647 


22 


1.1790267 


3.20 


2.49945 


- 0.606 


23 


1.1790301 


3.04 


2.49636 


- 0.564 


24 


1.1790316 


2.85 


2.49387 


- 0.518 


Expected 


1.1790302 


2.91 


2.49642 


-0.529 



Table 1. Values of /3 C , /3 c co//i, -B4 and Co obtained from linear fits with the 
(2(n-&)/3 ^ 2* 1 / 3 ) set of 6 intersections described in the text. 



We have also performed fits of the intersections for each value of N using all 
the possible intersections involving the 5 values of the linear size immediately below. 
This amounts to 15 intersections, just as for the determination of C . This procedure 
maximizes the number of data points for a given size interval and should be suitable 
when larger numerical errors are present. We have repeated the calculation with 
smaller and larger values of e. The results for the f3 c and B4 are shown in figure [U 
One sees that the results become erratic for e > 0.3. There is a line crossing that 
appears in both graphs for 0.1 < e < 0.15. We believe that this crossing reflects the 
compromise discussed above, even though it occurs at a value of e about two times 
larger than initially guessed. The estimates at the crossing are f3 c = 1.179025(5) and 
B 4 = 2.492(3). In general, we do not know precisely the crossing value and so larger 
error bars should be set by using the variation with e. From the figures, this should be 
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Model 


B 4 


Ref. 


Method 


Intersections 


Ising 3 


1.590 (15?) 


m 


MC + reweighting 


extrapolation 


Ising 3 


1.591 (9) 


® 


MCRG + reweighting 


(128,256) 


Ising 3 


1.604(1) 


® 


MC at best C 


- 


Ising 3 


1.603(1) 


m 


Improved action 




Ising 3 


1.603(1) 


A 


Improved action 




SU{2) N T = 4 


1.620(30) 




MC at various (3 


extrapolation 


SU{2) N T = 6 


1.570(80) 


55 


)! 


55 


SU(2) N T = 8 


1.520(100) 


11 


55 


55 


SU{2) N T = 2 


1.622(3) 


m 


MC at various (3 


(24,32) 


SU{2) N T = A 


1.602(27) 


ii 


55 


J) 


SU(2) N T = 8 


1.584(39) 


ii 


55 


)5 



Table 2. Values of B4 found in the literature for the two models discussed in 
the text. The last column provides additional information about a particular 
intersection (TV, TV') used or about the fact that an extrapolation has been made. 

at most 1 x 10~ 5 for j3 c and 0.005 for B4 for N — 256. The accurate values are clearly 
within the errors bar for the estimates obtained with the 6 and 15 points fits. It is 
also instructive to compare the 15 points fits in figure H] at e = 0.05 with the values 
of the 6 point fits given in table [T] for the same values of n. This gives an idea of the 
variability of the extrapolations as we change the data used to make the linear fits of 
equation (|2.8p . 

It is important to realize that near the value of e where the various curves of 
figure [4] cross, the extrapolations are approximately the same for all the values of n. 
This means that picking e properly allows to get optimally accurate estimates with 
volumes not too large. One should also appreciate that there are 35 values of e in each 
graph of figure |4] and that for each values of e, we need more than 100 calculations of 
B4 to determine the intersections. Doing the same number of calculations using the 
MC method for spin models or lattice gauge theory would take a very long time. 

5. Comments about the literature 

In this section, we comment on existing numerical results found in the literature. We 
first discuss existing estimates of B4 for the Ising universality class. The methods 
used can be divided into those that rely on intersections [H [6] , and those that attempt 
to work as closely as possible to the nontrivial fixed point [SJ [71 Note that [5] 
takes into account the nonlinear effect of /2 and also the effect of the second irrelevant 
direction. The results are summarized in table [2j The error bars in [3] are based on 
the visual estimate U = 0.47 ± 0.005 from figure 7. The agreement between the two 
results based on improved action that reduce the effect of the irrelevant variables [8] 
as well as the agreement with [5] indicates that 1.603 should be accurate with about 
one part in a 1000. The discrepancy with [5] seems to be of statistical origin: if we use 
the fit given in equation (14) of [5] together with N = 256 and \K - K c \ < 4 x 10~ 6 
we obtain an absolute error less than 10~ 3 for B4 if we neglect fi. 

We now discuss the nonlinear effects for specific calculations made with the fixed 
interval procedure. For the 3 dimensional Ising model, (3 C (usually denoted K c ) is 
known with great accuracy and consequently, most of the graphs showing B4 use a 
very narrow k interval where the differences between linear fits and numerical data are 
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Figure 4. Infinite volume extrapolations of /3 C and based on 15 point linear 
fits from the intersections among the B4 curves at N = 2"/ 3 and the 5 values of 
N immediately below, for n between 19 and 24. 
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Model 


Ising 3 


Villain XY 


SU{2) N T = A 


SU{2) N T = ± 




m 


m 


m 


EM 


Data used 


fig. 6, N = 16 


Fit on fig. 11 


fig. 9, N a = 8 


fig. 2, N a = 16 


hlh 


>9(3) 


3 


2 


> 3(1.5) 


\nN x l v \ max . 


5.5 


0.4 


0.47 


0.4 


^max. 


< 0.6(2) 


0.13 


0.24 


< 0.13(6) 



Table 3. Values of /1//2, \kN 1 / " \mtue. and E maJ . for the four references discussed 
in the text. 



difficult to see. For instance, if we use the N = 16 data in figure 6 of [2], it is possible 
to draw a line that stays within the error bars of all the data points. Consequently, we 
can only set a lower bound on /1//2 from this figure by assuming that Ao, defined in 
|Appcndix B| is less than the error bars that we estimated to be less than 0.02. This 
bound is consistent with the numerical estimate of [5]. 

A similar situation is encountered with figure 2 of [H] for SU(2) LGT with 
N T = 4. We focus our analysis on the N a = 16 data points. If we draw a line from the 
point with the lowest horizontal coordinate to the point with the largest one, the three 
points inside appear to be above the line. However, the line is within the three error 
bars. Using the data points kindly provided by the author, we obtained Ao ~ 0.01 
which is slightly less than the errors bars. The lower bound quoted in the table is 
based on Ao < 0.015. On the other hand, a clearly non-zero value for Ao can be seen 
in figure 11 of [4] and figure 9 of [10] and relatively stable values can be obtained for 
/1//2. If we now denote \nN x ^\ max , the maximum value of \kN 1 / v \ used in the linear 
fits made to determine the intersections (assuming that such a fits were performed, 
otherwise we rely on the width of the figure) in each reference, we can estimate the 
relative size of the nonlinear effects by calculating 

e m ax. = (f2/fl)\KN 1 ^\ max . . (5.1) 

The numerical values are shown in table [3l There are no large values of e max .- This 
is because when e max . ~ 1, the failure of linear fits is obvious. On the other hand, 
systematic errors due to unaccounted nonlinear effects may be reduced by using a 
smaller e max . with the procedure discussed in section [3l 

6. Conclusions 

We have proposed a new method designed to reduce possible nonlinear effects in the 
estimates of (3 C and B4. For the Ising hierarchical model with a volume corresponding 
to a linear size of 256 in 3 dimensions, the method gives results which agree with 
independent accurate estimates with a relative accuracy better than one part in 
100,000 for f3 c and 5 parts in 1000 for B4. The intrinsic nonlinear effects (measured 
in terms of /2//1) in this model are roughly of the same size as those found the other 
models that we discussed in section[S] Figures [5] and [3] shows the potentially disastrous 
effect of ignoring nonlinear effects. In graphs with large numerical errors, this effect 
may be overlooked. The small values of e max found in the literature indicate that 
nothing drastic should appear in the cases considered. However, it would be interesting 
to repeat and extend these calculations for larger volumes using the method proposed 
here. The optimal value of e can in principle be obtained by looking at crossings of 
extrapolated values. This suggests that the procedure should be followed for a few 
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different values of e. We are planning to consider the effects of statistical errors on the 
new method by applying it to the 3 dimensional Ising model and finite temperature 
gluodynamics. 
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Appendix A. Numerical Calculations 

For completeness, we explain how we calculated B4 for the hierarchical model and we 
give the details of the numerical calculations of the parameters. The interest of the 
hierarchical model is that we can blockspin exactly. The basic formula for the Fourier 
transform of the local measure is 

Rn+r(k) = C n+lC ^W^ + HoVo^ Rl{k) {A1) 

with 

c = 2 1 - 2 ' D . (A.2) 

Note that we do not rescale the field as in a renormalization group transformation. 
C„+i is a constant that we adjust so that i?„+i(0) = 1. When it is the case, we have 

00 

R n (k) =l + J2 (-k) 2q /(2q)l < (<Pn) 2q >n , (A.3) 
q=2 

with 4> n the sum of all the spins. If we write 

R n {k) = 1 + a 7hl k 2 + a ni2 k 4 + ... + a n .i malc k 21 *™* . (A.4) 

Then 

B 4 = Qa n ^/a 2 n ^ . (A. 5) 

The critical values are obtained by plugging the values of a\ and from the nontrivial 
fixed point. See [52] for details and references. The value of Co +ci{kN 1 /' j ) is obtained 
by increasing N but keeping kN 1 ^ constant. The asymptotic values for different 
kN 1 ^ can be fitted very well with a line which gives Co and c\. We then subtracted 
the effect of c + c^nN 1 /") and obtained /1 and ji using discrete derivatives near 
K = 0. 

We used calculations at very large volume (N ~ 10 6 ). In the text, we show that 
these accurate results can be reproduced with a new procedure using data at much 
smaller volume. From Refs. [2171 121) . 

v = 0.649570365 , 

uj = 0.655736 , 

Bi(j3 c , 00) = 2.49641845 , (A.6) 

P c = 1.17903017044626973251 . 

Using these values, and numerical results a fixed values of kN 1 ^ and large TV, we first 
determine the terms and find 

c = - 0.529 , 

ci = - 0.236 . (A.7) 
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Subtracting these effect and taking discrete approximation of the derivative with 
respect to k near 0, we find 



One should appreciate that we have been able to determine numerically the 8 
parameters of equation (|2.7| with good numerical stability. 

Appendix B. Estimate of the nonlinear effects with pencil and ruler 

In this appendix, we explain how to estimate the order of magnitude of /1//2 from 
a graph illustrating approximate data collapse such as figure [1] using a pencil and a 
ruler. A subset of the data of figure [1] is shown in figure iBll together with lines that 
can be drawn on the original graph. The line joins two points with approximately 
opposite x coordinates. We denote these two points (x,yi) and (— x, 1/2)- We call Ao 
the difference between the line and the data (assumed to be a quadratic function) at 
. It is then easy to show that 



In figure [BT1 we have approximately x ~ 1.24, y^ — 2.66, yi ~ 2.06 and Aq ~ 0.07, 
which gives the estimates fi ~ —0.24 and — —0.045 which are in reasonably good 
agreement with the previous estimates. Note that we are mostly interested in the ratio 
I/1//2I and that this ratio is invariant under a multiplicative rcscaling. This is useful 
if we are considering graphs with other quantities displayed (for instance 1 — Bi/3). 
For the model considered here, I/1//2I — 5. 

The above estimates are based on the data with N = 32. In figure [BlT we also 
see part of the data for N = 64 and we see that it is not much above the N — 32 data. 
From equation (|2.7[) . it is clear that if we neglect c\, the entire curve is translated 
uniformly in the vertical direction. But the estimates of fx and /2 are based on 
differences and do not depend on this translation provided that we only use data for 
one value of N. This is confirmed by using the N — 16 data as shown in figure IbTI 
where we obtain an estimate of /1//2 very close to the one quoted above for N = 32. 
This shows that up to a subtraction and rescaling of the horizontal coordinate, we 
can rely on a graph of B4 versus (3 at relatively small volume. It is clear that the 
subtraction of (3 C requires to have an estimate of this quantity. However, this estimate 
does not need to be very precise to get the order of magnitude of /i//2- 
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